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1. Introduction 

The study of type la supernovae led to the 1998 discovery of the acceleration of 
the cosmic expansion [TJ, which has prompted theoretical physicists to look for an 
explanation. Many models of this accelerated dynamics postulate the existence of dark 
energy, a diffuse mysterious form of energy with exotic equation of state P ~ — p which 
amounts to 70% of the critical energy density, with the remaining 30% comprised of 
dark and ordinary matter. As an alternative to exotic dark energy models, it has 
been proposed that infrared modifications of gravity at the largest scales could be the 
explanation of the cosmic acceleration [2]. The idea is either to consider the Einstein- 
Hilbert actional with cosmological constant 

S EH = [ d A x^g~(R - 2A) + S {matter) (1.1) 



2k „ 

or some form of dark energy mimicking A, or to change gravity as described by the 
alternative action 

S=^J d 4 xy/=g~f(R) + S {matter) , (1.2) 

where f(R) is a non-linear function of R. This "f(R) n or "modified" gravity has a 
long history (see Ref. [4] for an historical perspective): it was originally conceived 
as a mathematical alternative to Einstein's theory [5j; quadratic corrections to the 
Einstein-Hilbert action were discovered to be necessary for the renormalization of general 
relativity [6], and then used to generate inflation in the early universe without scalar 
fields [7\; see Ref. |8| for a more modern perspective. 

Modified gravity comes in three possible forms: the first is metric f(R) gravity, in 
which the metric is the only independent variable and the variation of the action (11.21) 
with respect to g ab yields the fourth-order field equations 

f(R)R ab - ^-g ab = V a V b f{R) - g ab nf(R) + kT™ . (1.3) 

The second possibility, named "Palatini f(R) gravity" consists of treating the metric 
and the connection as independent variables (i.e., the connection T^ c is not the metric 
connection {£ c } of g ab ). If the (non-metric) connection is assumed to be symmetric, the 
field equations (which are now of second order) are [9] 



f{R 



f[R)R ab - ^-^g ab = nT^ , (1.4) 



-gf{R)g ab +V d ^- g f(R)g^ 5? = 0, (1.5) 



where R ab is the Ricci tensor of the non-metric connection, R = g ab R ab , and V c is the 
covariant derivative operator of T^ c . 

| Here K = 87rG, R is the Ricci curvature, A is the cosmological constant, and g is the determinant of 
the metric tensor g ab . We follow the notations of Ref. [3]. 
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The third possibility, "metric-affine gravity", allows also the matter action to depend 
on the independent connection, resulting in more general field equations |10| . In this 
paper, we only consider the metric and Palatini formalisms. 

First, let us focus on metric f(R) gravity. In order for these theories to be viable, 
several criteria must be met, which have been analysed in the recent past. The theory 
must have the correct Newtonian and post-Newtonian limits [HI [12] and it must not have 
short timescale instabilities, which is achieved if /" > |T3l [EH [EU \W\ (see Refs. [17] 
for studies of various instabilities). The theory must possess the correct cosmological 
dynamics, i.e., an inflationary era, followed by a radiation era and then a matter era 
and, finally, an acceleration era. Furthermore, there must be smooth transitions between 
consecutive eras, which is not always achieved, with many models being ruled out [18j. 
The theory must also be ghost-free, which is the case for f(R) and for Gauss-Bonnet 
gravity, but not for theories of the form / (R, R a bR ab , R a b ca 'R abcd ) pEH]. Another criterion 
for validity that should certainly be satisfied is that the theory possesses a well-posed 
initial value formulation. Here we focus on that aspect. In addition to being a viability 
criterion, the Cauchy problem proves very useful in the numerical integration of the field 
equations. 

The initial value problem was studied for theories of the form 

S = J d 4 x^9 (j^ + aR ab R ab + (3R 2 + jR^B^j + S matter (1.6) 

in Refs. [20, EE] with the conclusion that the Cauchy problem is well-posed. Due to 
the Gauss-Bonnet identity in four spacetime dimensions, the Kretschmann scalar can 
be dropped from the action. In the following, we consider the Cauchy problem for the 
theory described by the action (11.21) and for a general form of the function f(R), without 
restricting to the specific theory f(R) — R + aR 2 considered in Refs. [2Ql EE]- Instead 
of proceeding directly on the field equations (11.31) . we will use a well known dynamical 
equivalence relation mapping f(R) into scalar-tensor theories [22], EE] and we will address 
the Cauchy problem in scalar-tensor gravity, relying on recent results on this subject |23j. 
Our main result is that the Cauchy problem for metric f(R) gravity is well- formulated, 
and is well-posed in vacuum. For Palatini f(R) gravity, instead, the Cauchy problem 
is not well-formulated nor well-posed due to the presence of higher derivatives of the 
matter fields in the field equations and to the impossibility of eliminating them. Here, 
the system of 3 + 1 equations of motion is said to be well-formulated if it can be re- written 
as a system of only first-order equations (in time and space) in the scalar field variables. 
When this set can be put in the full first order form 

d t u + M i V i u = S(u) , (1.7) 

where u collectively denotes the fundamental variables hij,Kij, etc. introduced below, 
M % is called the characteristic matrix of the system, and S (u) describes source terms 
and contains only the fundamental variables but not their derivatives. Then, the initial 
value formulation is well-posed if the system of partial differential equations is symmetric 
hyperbolic (i.e., M l are symmetric) and strongly hyperbolic (i.e., SiM 1 has a real set 
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of eigenvalues and a complete set of eigenvectors for any 1-form Sj, and obeys some 
boundedness conditions). We refer the reader to Ref. [24] for more precise definitions 
and for a complete discussion. 

The plan of this paper is the following: in Sec. 2 we recall the equivalence between 
f(R) and scalar-tensor gravity. In Sec. 3 we proceed to show that the Cauchy problem 
is well- formulated for the scalar-tensor equivalent of metric f(R) gravity. The scalar- 
tensor equivalent of Palatini f(R) gravity is considered in Sec. 4, and Sec. 5 contains a 
discussion and the conclusions. 



2. The equivalence between f(R) and scalar-tensor gravity 

We briefly review the equivalence between f(R) and scalar-tensor gravity, starting with 
the metric formalism. By introducing the extra scalar field 0, the modified gravity action 

S=^J d 4 x^f(R) + S (matter) (2.1) 

can be rewritten as [22], [21] 

S = / d 4 x^g~ [tp((p)R - V((j))] + S {matter) (2.2) 
2k J 

when f"(R) 7^ 0, where 

i/> = f(0) , V{<j>) = 0f (0) - f{<j>) . (2.3) 

The action (12.21 ) coincides with ( 12.11) if = R. Vice-versa, the variation of ( 12 .21) with 
respect to yields 

"%-% = i"-+)rw = °, (2.4) 



which has no dynamical content but implies^ that = R when f"(R) ^ 0. (12.21) is a 
Brans-Dicke action of the form |26| 



S=^J d 4 xy/=g~ ipR - |V C ^V C ^ - Uty) + S (matter) (2.5) 

with Brans-Dicke field ip, U(tp) = V[(p(tp)], and Brans-Dicke parameter lu — 0. This 
is called O'Hanlon theory or "massive dilaton gravity" and was originally introduced in 
order to produce a Yukawa term in the Newtonian potential when taking the Newtonian 
limit [27J. The field equations are 

G ab = - T { ™ ] - -L U^)g ab + - ( VaV 6 V - <M&A) , (2.6) 

3Q0 + 2U (tp) — tp = k . (2.7) 

dtp 

Let us consider now modified gravity in the Palatini formalism. The Palatini action 

S=^J d 4 x^f(R) + S {matter) (2.8) 

§ An action of the form (|2.ip with the property that f"(R) ^ is sometimes called "R- regular" [25] . 
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is equivalent to 

S 



1 

2k 



d 4 Xy/—g 



fix) + fix) (R-x)]+S^ atter \ 



(2.9) 



the variation of which with respect to \ yields \ — R- One now uses the field (ft = fix) 
and takes advantage of the fact that the curvature R can be seen as the Ricci curvature 
associated to the new metric 



lab 



f'(R)g ab (2.10) 

conformal to g a b [2H [10]. The two conformal frames are physically equivalent (see the 
recent extensive discussion of |29| and references therein) . The transformation property 
of the Ricci scalar under conformal transformations is [30], [3] 



R = R+— V c 0V, 

2(f) 



3 -a, 

2 v 



which allows the action 02.9jl to be rewritten as 



S ~ 2k 



d A X\J—g 



+ ^VWc0-^(0) 



| (^{matter) 



(2.11) 



(2.12) 



where an irrelevant boundary term has been omitted and 

V(^)=M<P)-fix(<P)] ■ (2-13) 

The action (12 . 121) describes a Brans-Dicke theory with parameter uj = —3/2. This theory 
is generally regarded as a pathological case [HI], but is sometimes studied [32] . 
A general Brans-Dicke theory 



S 



BD 



1 
2~K 



d Xa 
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<j)R - -V c 0V c - V((f)) 



I ^(matter) 



(2.14) 



leads to an effective gravitational coupling appearing in a Cavendish experiment (and 
in the theory of cosmological perturbations [33] ) 

2 (u + 2) 1 



G 



eff 



(2lu + 3) <p ' 

which becomes ill defined for u = —3/2. Similarly, the equation of motion for 



(2lu + 3) Op = K T (m) + 0^ - 2V(4>) 



(2.15) 



(2.16) 



reduces to an algebraic identity when uj = — 3/2 and the field becomes non-dynamical. 
This has obviously consequences for the initial value problem of the theory, which are 
discussed later. 

Our strategy will be to study the Cauchy problem of f(R) gravity by reducing it 
to the initial value problem of lj = 0, —3/2 Brans-Dicke theory. Nevertheless, we start 
by considering more general scalar-tensor theories in the next section. 
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3. The initial value problem of scalar-tensor gravity 



The Cauchy problem of Brans-Dicke and scalar-tensor gravity theories was considered 
in Refs. [20l EH, [34]. Noakes [20] showed that the Cauchy problem is well-posed for a 
non-minimally coupled scalar field with action 



S 



d A x\J—g 



1 

2k 



i2__v^v c ^-y(0) 



(3.1) 



without matter. Originally, this theory was introduced by regarding as a form of 
(quantum) matter on a fixed curved background, not as a gravitational scalar field akin 
to the Brans-Dicke field, nor as a source of gravity on the right hand side of the Einstein 
equations. There is little doubt that this is the reason why Noakes did not include 
"ordinary" matter in the theory, which would complicate the discussion and weaken his 
conclusion of well-posedness of the Cauchy problem [JJ. However, ( 13.11) can be legitimately 
regarded scalar-tensor action. 

Cocke and Cohen [34] sketched a study using Gaussian normal coordinates for 
Brans-Dicke theory with a free Brans-Dicke scalar. A covariant and systematic approach 
to the Cauchy problem of scalar-tensor theories of the form 



S 



d A Xy 



-g 



MR 

2k 



-V c 0V c 0-V(0) 



I ^(matter) 



(3.2) 



has recently been advanced by Salgado [23], who showed that the initial value problem 
for these theories is well-posed in the absence of matter and well-formulated otherwise. 
For our purposes it is necessary to generalize Salgado's results to more general scalar- 
tensor theories of the form 



s 



d Xy 



-g 



- V c 0V c - V{(j>) 



I ^{matter) 



(3.3) 



2k 2 

which contain the additional coupling function u) (</)). In the following, we show in 
detail how this generalization proceeds and then specialize the discussion to the cases 
uj = 0,-3/2 which are the main focus of our interest. We follow as closely as possible 
the notations of Ref. [23] in order to facilitate comparison. 

The field equations are (we set k = 1 in the rest of this paper) 
1 



G n 



1 

+ 7 



/ 



[/" (V a 0V 6 - g ab V c <pV c <P) + f (V a V b <j) - g ab U^ 



^ V a 0V fe - -£ afe V c 0V c - V{4>)g ah + T L 







" ab 



^□0+^i?-F'(0) + yV c 0V c 



(3.4) 



(3.5) 



|| This example shows that there is no clear-cut distinction between "gravitational" and "non- 
gravitational" scalar fields in relativistic theories of gravity (see Ref. |35J for a discussion and 
implications). 
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where a prime denotes differentiation with respect to 0. Eq. (13.41) is rewritten as the 
effective Einstein equation 

T (eff) 



G n 



where the effective stress-energy tensor is [23] 



<eff) 



rp{f) I rp(4>) | rp(m) 

1 ab ' 1 ab ' 1 ab 



(3.6) 
(3.7) 



with 

T% } = f"(4>) (V a 0V b - ^V c 0V c 0) + f{4>) (V a V fe - g ab D4>) , (3.8) 
which contains second order derivatives of <t> and is identical to the one in Ref. 1 23 1 , and 



= w(0) ( V a 0V 



-^V c </»V c - V{(t))g ab 



(3.9) 



which has canonical structure (i.e., quadratic in the first-order derivatives of <fi with no 
second order derivatives) and differs from the one of Ref. [23] by the presence of the 
coupling function cu(0). By taking the trace of eq. (13.61) and solving for □</>, one obtains 



□c 



^ 2fV(<P)+fV'( ( f>) + 


"-^-4(^ + 3/")" 




f 


[., , 3(/') 2 " 





.(3.10) 



The 3 + 1 ADM formulation of the theory proceeds by defining the usual quantities lapse, 
shift, extrinsic curvature, and gradients of (3j EH [23]. Assume that a time function t 
exists such that the spacetime (M,g ab ) can be foliated by a family of hypersurfaces E t 
of constant t with unit timelike normal n a . The 3-metric is defined by h ab = g ab + n a n b 
and h a c is the projection operator on E t , with 

/i a6 n b = /i ab n a = , hah bf , = h n 



(3.11) 
(3.12) 
(3.13) 
(3.14) 
(3.15) 

- bl . . . n-iv f -j. - di... (3.16) 
for any 3-tensor ( 3 'T ai ~ bl .., , with Dih ab = 0. The spatial gradient of the scalar field is 



n U a = — 1 , l lab' l = HabU = u , lb a IL bc 

The metric is written using the lapse function N and the shift vector N a as 
ds 2 = — (iV 2 — N { Ni) dt 2 - 2N i dtdx i + h ij dx i dx j 
= 1,2,3), with N > 0, n a = -NVJ and 
N a = -h a b t b , 

where t a is a time flow vector satisfying t a V a t = 1 and 

t a = _ N a + Nn a 

so that N = —n a t a and N a n a = 0. The extrinsic curvature of T, t is 

K ab = -h a c h b d V c n d 
and the 3D covariant derivative of h ab on S t is defined by 



D (3) T a 



h ai Cl . . . h d \, . . . h s M! ] T c 



Qa = D a 



(3.17) 
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while its momentum is 

n = £ n = n c V c (3.18) 



with 



Kl3 = -V inj = -— -Ji + DM + DjNi , (3.19) 



1 / dhij 
2N \~dt 

U = ^(d t <p + N C Q C ) , (3.20) 



and 



dtQi + N l diQi + QAN 1 = A (NIL) . (3.21) 
The 3 + 1 decomposition of the effective stress-energy tensor follows: 

T a ( fe e//) = j (S ab + J a n b + J b n a + £n n n 6 ) , (3.22) 



where 



a -l cl drp{eff) _ 1 / rr(/) Q 0) qM"\ /„ O o\ 

>3a6 = "a Kb 1-cd -J y^ab + ^afe + ^ab j i (d.2dj 

Ja = ~h a c T^n d = i (JW + JW + JH) , (3.24) 

£ = n a n b T^ ff) = j (E {f) + # + £ (m) ) , (3.25) 

and T< e ^ = S — E, where = T( e ^)° a and 5 = S a a . By means of the 

Gauss-Codacci equations, the effective Einstein equations projected tangentially and 
orthogonally to T, t provide the Hamiltonian constraint [HI [23] 

(3) i? + K 2 - KijK* = 2E , (3.26) 
the vector (or momentum) constraint 

DtK l i - D t K = Ji , (3.27) 
and the dynamical equations 

dtlCj + N l diK l j + K\djN l - K l j diN i + DWjN 

RjN - NKK'j = — [(S- E) 5) - 2Sj] , (3.28) 

where K = K\. The trace of this equation yields 

d t K + N l diK + (3) AiV - NKijK ij = — (S + E) , (3.29) 

where ® A = D l Di. As remarked in [20], [23], the presence of second order derivatives 
of the scalar (seen here as an effective form of "matter") could potentially render the 
Cauchy problem ill- formulated. However, this is not the case because these derivatives 
can be eliminated in most cases (with one notable exception explained later) |23| . 



The Cauchy problem of f(R) gravity 9 
The /- and 0-quantities appearing in eqs. (13.23l) -( l3T25l ) are then calculated, yielding 



E^ = f (D C Q C + Kit) + f"Q 2 , 
JW = -f (KQ C + DJl) - f"UQ a 



(3.30) 
(3.31) 



S$> = f (D a Q b + UK ab - h ab U<t>) - f" [h ab (Q 2 - IT 2 ) - Q a Q b ] , (3.32) 



where Q 2 = Q C Q C . Other useful quantities are 

£(/) = /' (D C Q C + KU- 30(f)) + f" (3U 2 - 2Q 2 ) 
S (f) - #(/) = -3 f'U(j> - 3 /" (Q 2 - n 2 ) . 



(3.33) 
(3.34) 



All these /-quantities coincide with those of Ref. |23j. However, the following 0- 
quantities differ from the corresponding ones of [23] by the presence of terms in ui and 



eP> = | (n 2 + q 2 ) + v{<j>) 



= ooQ a Q b - h 



ah 



UJ 



(q 2 - n 2 ) + v{4>) 



while 



and 



(3.35) 
(3.36) 
(3.37) 

(3.38) 



The "total" quantities entering the right hand side of the 3 + 1 field equations are then 

E = j[f ( DC Qc + KU ) + (/" + 1) Q 2 + | n2 + v (<t>) + E(m) ] > ( 3 - 4 °) 



J a = j [-/' (K a c Q c + D a Tl) - (/" + lu) TlQ a + 4™)] 



(3.41) 



S ab = j {/' (D a Q b + UK ab ) - h ab [(/" + |) (Q 2 - n 2 ) + K(0) + f □</>] } 



1 

7 



^ + /") Q a Q b + S, 



(m) 



(3.42) 



while 



s = jlf ( DC Qc + UK ) ~ W (<P) - 3 fa, 



i 

7 



^ -2 Q (fit , W \ n 2 _j_ ^(m) 



(3.43) 



S-E = 
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\ [(3f + u) (n 2 - Q 2 ) - 4V(4>) - 3f'D(f> + S {m) - £ (m) ] , (3.44) 



/ 



S + E = - [2f (D C Q C + KTL) - f'Q 2 + (3/" + 2w) n 2 ] 
+ j (-2V{(j>) - 3f'D(f> + S (m) + £ (m) ) . 



(3.45) 



The Hamiltonian constraint becomes 

{3) R + K 2 - KaK ij - - 

f 



f (D C Q C + KU) + |n 2 + ^ (u + 2f") 



= - (E™ + V (<!>)) , 
the momentum constraint ( 13.271) is 



(3.46) 



DiK l i -D t K+- [f {Ki c Qc + An) + (u + /") ng,] 



J 



(m) 



the dynamical equation (13.281) is written as 

dt&j + tffylPj + K\djN l - K/diN* + D'DjN - (3) R^N - NKICj 



(3.47) 



+ Yf t f " ( g2 ~ n2 ) + 2V{(f)) + fn ^ % + T" + UKij ^ 

+ ^r(u + f") Q l Q 3 = ^ \(S (m) - 5) - 2S {m) y , 
J 2/ 



(3.48) 



with trace 



d t K + N l diK + (3) AN - NK^K^ - ^-j- (D C Q C + UK) 

+ ^ [f'Q 2 - ( 2uj + 3 f") n2 ] = Yf (~ 2V ^ ~ 3/ ' D0 + s(m) + E{r t¥ 9 ) 

where (cf. Ref. [23]) 

CJl - nAT - Q C D C (lnJV) - £> C Q C = -□</> 



3(/Q 2 
2/ 



/ 



3(/T 
2/ 



2 



- 2/V(0) + fV'{4>) 

r 



--^ + 3/0- 



V c 0V c 



(3.50) 



In vacuo, the initial data (hij, Kij, (j), Qi, IT) on an initial hypersurface So must satisfy 
the constraints (13.461 ) and ( 13.471) plus 

Qi - Di<f> = , (3.51) 
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(3.52) 



In the presence of matter, the variables E( m \ ji m \ must also be specified on the 
initial hypersurface. Gauge-fixing is equivalent to assigning lapse iV and shift N a (various 
gauge conditions used in numerical analysis are reviewed in [23]). The system (13.460 - 
( 13.491) contains only first-order derivatives in both space and time once the d'Alembertian 
□0 is written in terms of 0, V c 0V c 0, /, and its derivatives by using eq. ( 13.501) . This 



system differs from the corresponding one in |23| only by terms in oj((p) and its first 
derivative. From now on, everything will proceed as in Ref. [23] (with one exception 
discussed in the next section). The reduction to a first-order system indicates that the 
Cauchy problem is well-posed in vacuo, this time also for the more general scalar-tensor 
theories containing the extra coupling function w(0), and that it is well- formulated in 
the presence of matter. We do not repeat Salgado's analysis here, referring the reader 
to [23]. We are now ready to translate the results in terms of f(R) gravity. 



4. The Cauchy problem of modified gravity 

In the notations of Ref. [23] that we adopted, Brans-Dicke theory [26J is described by 
u;(0) = ou /(p with u the constant Brans-Dicke parameter, /(0) = 0, and V — > 2V. 
This yields the constraints 



{3) R + K 2 - K {j K 



D c Qc + KU + ^ (II 2 + Q 2 ) 



- + V 



(4.1) 



DiK l i — DiK + ^ \^K i l Q l + An + ^nQiJ 



J, 



(m) 



and the dynamical equations 

d t K l j + N l diK\j + K\djN l - K/diN* + D^jN 



(4.2) 



- (3) R i jN - NKK'j + ^-8) (2V(4>) + D<f>) + j (D'Qj + UK^) 

+ ^Q l Qi = ^- ((s {m) - £ (m) ) 5) - 2s^ y , 



(4.3) 



with 



d t K + N l diK +^ AN - NKijK ij - ^ (D C Q C + UK) - ^IT 2 

<p <p A 



[-2V(<f>) - 3D0 + S {m) + E (m) ] 



N 
20 



2V(0) + 0V(0) + ^(n 2 -g 2 ) . 



(4.4) 



(4.5) 
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>From the discussion of the previous sections, we conclude that metric f(R) gravity, 
which is equivalent to ujq = Brans-Dicke gravity, has a well-formulated Cauchy problem 
in general and is also well-posed in vacuo. On the other hand, Palatini f{R) gravity, 
which is equivalent to a Brans-Dicke theory with ojq = —3/2, is in trouble. In fact, for 
this value of the Brans-Dicke parameter, the d'Alembertian □</> disappears from eq. (14.51 ) 
and the field is not dynamical — it can be arbitrarily assigned on a region or on the 
entire spacetime, provided its gradient satisfies the degenerate equation (14.51) . which 
reduces to a constraint. 

As a consequence, in the Palatini version, it is impossible to eliminate □</> from 
the system (14. II )- (14 Afi unless D0 = 0. This condition includes both the possibility = 
constant, in which case the theory degenerates to general relativity, and the case of a 
harmonic 0, which is associated to a null 0-wave. Apart from these special cases, the 
previous considerations constitute a no-go theorem for Palatini f(R) gravity, which has 
an ill-formulated Cauchy problem in vacuo and, therefore, can hardly be regarded as 
a viable theory. The absence of a well-posed Cauchy problem was briefly noticed in 
Ref. [20]. 

5. Discussion and conclusions 

While metric f(R) gravity has a well- formulated initial value problem (and a well- 
posed one in vacuo), Palatini modified gravity does not. In conjunction with recent 
criticism of the Palatini formalism (according to which the correct Palatini variation 
necessitates a Lagrange multiplier method which is usually neglected and, when taken 
into account, produces essentially the same theory as the metric version [38] — see also 
[39]). this fact probably leads to the demise of Palatini f(R) gravity as a legitimate 
model of the cosmic acceleration. However, something can be learned from the situation 
of the Palatini version of modified gravity. It is interesting to speculate on the possible 
physical consequences of the failure of a theory to have a well-posed initial value problem 
for the reasons discussed in the previous section. As seen above, the difficulty originates 
from the fact that the effective "matter", i.e., the scalar field 0, is described by a 
stress-energy tensor containing derivatives of of second order instead of first. This 
non-canonical form makes it very difficult, if not impossible, to satisfy the energy 
conditions [3]. This feature in itself should not necessarily be regarded as a curse or 
failure. However, while for the lu = equivalent of metric f(R) gravity the second 
derivatives contained in D0 can be eliminated, for the uj = —3/2 equivalent of Palatini 
f(R) gravity this is not possible and, when integrating twice to solve the field equations, 
the solution for the metric will depend on 0, not only on an integral of (for example 
through the usual Green function integral). Then, the metric will be very sensitive to 
variations in that are not smoothed out, or averaged, by an integral. This extreme 
sensitivity, contrary to ordinary situations, may spoil the continuous dependence on 
initial conditions. Physically, this is the situation encountered in modelling stars in 
Palatini f{R) gravity. In Ref. [37] it is found that the dependence of the metric on 
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higher order derivatives of the matter fields makes it so sensitive to small variations that 
it is impossible to build even a model of a polytropic star, while a polytropic gas is a 
perfectly reasonable form of matter even in a Newtonian star. It is auspicable that this 
dependence on higher order derivatives of the matter fields be avoided in future theories 
of gravity, for example, by including higher order terms in R a bR ab in the action. 
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